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Effect  of  Polydispersity  on  the  Cloud  Point  Curves 
of  Polymer  Mixtures 

Ryong-Joon  Roe  and  Learong  Lu 

Department  of  Materials  Science  and  Metallurgical  Engineering 
University  of  Cincinnati,  Cincinnati,  OH  H5221-0012 

AfimACX 

A  method  is  presented  for  the  calculation  of  cloud  point 
curves  of  polymer-polymer  mixtures  when  the  polymers  involved  are 
polydisperse.  The  method  is  based  on  the  Flory-Huggins  free 
energy  of  mixing  with  a  concentration-independent  X  parameter. 
Numerical  results  are  given  for  the  cases  where  the  molecular 
weight  distributions  are  represented  by  the  Schulz-Flory  type. 
When  the  two  polymers  have  similar  average  molecular  weights  and 
polydispersities,  the  cloud  point  curves  become  flatter  as  the 
polydispersity  increases.  When  the  two  polymers  have  similar 
average  molecular  weights  but  differ  in  their  polydispersities, 
the  cloud  point  curves  become  more  skewed  as  the  difference  in 
the  polydispersity  increases.  The  results  point  out  that,  if  the 
polydispersity  effect  is  not  properly  accounted  for,  the  value  of 
X  deduced  from  experimental  cloud  points  is  liable  to  be  in 
error,  especially  with  regard  to  its  temperature  coefficient  or 
its  concentration  dependence. 
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As  Is  well  known,  the  solubility  between  two  polymers  is  in 
general  severely  limited  because  of  the  rather  small  entropic 
gain  achieved  on  mixing  large  molecules.  Since  the  molecular 
weight  heterogeneity  modifies  only  the  entropic  term,  which  is 
already  small  to  begin  with,  it  could  be  expected  that  the 
polydispersity  influences  the  miscibility  behavior  of  polymers 
only  weakly.  This  is  certainly  true  when  the  polymers  involved 
are  all  of  extremely  high  molecular  weights.  In  many  recent 

1  o 

studies  J  of  thermodynamics  of  polymer  mixtures,  especially  of 
those  exhibiting  UCST  behavior,  however,  the  polymers  involved 
were  of  relatively  low  molecular  weights,  so  that  a  proper 
consideration  would  have  to  be  given  to  the  possible  effect  of 
polydispersity.  With  systems  exhibiting  LCST  behavior  the  phase 
separation  is  brought  about  usually  by  a  change  in  the  sign  of 
the  x  parameter  accompanying  a  temperature  change.  In  such  cases 
the  absolute  magnitude  of  the  X  term  is  small  and  therefore  the 
modification  of  the  small  entropic  term  by  the  polydispersity 
effect,  even  when  the  molecular  weights  are  high,  may  appreciably 
affect  the  observed  cloud  point  curves. 

Extensive  analyses  have  been  made,  by  Koningsveld,1*  Sole,® 
and  others,®  on  the  effect  of  polydispersity  on  phase  separation 
of  polymer  solutions,  i.e.,  solutions  of  a  polydisperse  polymer 
in  a  low  molecular  weight  solvent.  These  analyses  all  start  from 
the  Flory-Huggins  free  energy  of  mixing  (or  its  modification). 

In  this  work  we  likewise  take  the  Flory-Huggins  expression  for 
the  mixing  of  two  polymers  as  the  starting  point,  and  investigate 


numerically  the  expected  binodal  curves  as  a  function  of  the 
polydlspersity  of  the  two  polymers.  As  has  been  stressed  in 
recent  years, '  the  mean-field  approximation  embodied  in  the 
Flory-Huggins  expression  is  fairly  satisfactory  in  describing  the 
thermodynamics  of  polymer-polymer  mixtures.  A  recent  study  by 

D 

Joanny0  of  the  critical  behavior  of  polymer  mixtures  explicitly 
confirms  this.  This  is  to  be  contrasted  with  the  fact  that  for  a 
polymer  solution  the  Flory-Huggins  formula  becomes  grossly 
inadequate  near  the  critical  point  and  at  dilute  concentrations. 

Both  for  polymer-solvent  and  polymer-polymer  systems  the 
effect  of  polydlspersity  on  the  spinodal  curve  and  the  critical 
point  can  be  deduced  fairly  easily.  When  the  x  parameter  is 
independent  of  molecular  weights  (but  may  be  dependent  on  the 
concentration),  it  has  been  known^  that  the  spinodal  curve 
depends  only  on  the  weight-average  molecular  weight  and  the 
critical  point  on  the  weight-average  and  z-average  molecular 
weights.  No  such  simple  results,  however,  arise  in  the  case  of 
the  binodal  curve,  which  depends  on  the  details  of  the  molecular 
weight  distribution  and  can  therefore  be  deduced  only 
numerically. 

MEIliflJB 

We  start  the  analysis  with  the  Flory-Huggins  free  energy  of 
mixing  in  the  following  form: 

AGfl/ RT  =  £ ($ ip^P^ln^ip  +  2 +  i  (1) 

where  is  the  volume  fraction  of  the  p-mer  of  polymer  1, 
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and  x  is  a  function  of  temperature  but  is  independent  of  4>.  The 
chemical  potential  is  derived  from  equation  (1)  as: 

Aulp/RT  S  ln<|>1p  +  1  -  P(  <,i/Pin+<,>2/p2n)  +)(P^2)2  (3) 


where  p^n  is  the  number-average  DP  of  polymer  1.  When  two  phases 
coexist,  the  chemical  potentials  of  each  species  in  the  phases  A 
and  B  are  equal:  Ay^p  =  Ay®p,  and  it  therefore  follows  that: 

*V<?P  =  exp(olp);  *2p/(p2p  =  exP(°2p)  (i° 
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At  the  cloud  point,  at  which  an  infinitesimal  volume  of 
phase  A  has  just  separated  out,  the  composition  of  phase  B  is 
essentially  the  same  as  the  composition  of  the  whole  mixture: 


B 

^  1  p  p  (cloud  point)  (7) 

Eliminating  X  from  (5)  and  (6)  and  using  the  condition  (7),  we 
find: 

<*1*VV*  ♦  (*2*VV2  *  (♦l/ptn**2/p2n)  '  (<l’l/p1r,*Vp2n)(e) 


If  w<|(p)  and  w2(p)  are  the  (normalized)  molecular  weight 


*  % 
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distributions  (by  weight)  of  polymer  1  and  polymer  2 
respectively,  the  quantities  in  equation  (8)  pertaining  to  the 
incipient  phase  A  can  be  expressed  as: 

<£ |  =  4>-|  £w.|(p)  expfojp)  (9) 

P 

<f>^/p*n  =  <p >|  £  (7/p)  w.,(p)  exp(a1p)  (10) 

P 

and  similarly  for  ^  and  (f^/p**  Thus,  equation  (8)  can  be 
regarded  as  an  implicit  equation  fta^Og)  =  0  for  two  unknowns 
and  o 2 .  Another  implicit  equation  for  and  arises  from  the 
condition  that: 

<f>*  +  <f>£  =  1  (ID 

in  conjunction  with  equation  (9)*  Thus,  equations  (8)  and  (11) 
constitute  a  set  of  two  simultaneous  equations  which  can  be 
solved  for  and  for  given  values  of  4>i,  w-j(p),  and  w2(p).  X 
can  then  be  evaluated  from  either  equation  (5)  or  (6). 

For  numerical  calculation  of  the  cloud  point  curve,  we 
assume  the  molecular  weight  distribution  w(p)  to  be  given  by  the 
so-called  Schulz-Flory  distribution: 

w(p)  =  (ak+1pk/kl)  exp  (-ap)  (12) 

where  the  parameters  a  and  k  are  related  to  the  number  average 
DP,  Pn,  and  the  polydispersity  index,  8  =  Pw/pn,  by 

Pn  s  k/a;  g=i+  1/k  (13) 

The  summations  indicated  by  equations  (9)  and  (10)  can  then  be 
easily  evaluated  by  Integration.  The  numerical  solution  of  the 
equations  (8)  and  (11)  for  and  o2  is  then  straightforward. 


All  results  shown  below  were  obtained  with  a  PDP  11/23  computer 


RISULIS  ME  DISCUSSION 

Fig.  1  shows  the  cloud  point  curves  calculated  for  the  case 
where  plw  =  p2w  =  1000  and  B1  =  62  =  1.0,  1.1,  and  2.0 
respectively.  The  three  curves  go  through  the  common  critical 
point  at  =  0.5.  When  x  is  independent  of  composition,  the 
spinodal  curve  is  given  by: 

1/p1w<h  +  1/p2w*2  =  2X  (1«) 

and  the  critical  point  by: 

1/f2c  =  1  *  P2kP,VpiuP2z  <'5) 

and 

2XC  =  (1/p^  +  1/p^Kp^/p^  +  P2z/p2w)»  (16) 

from  which  it  can  easily  be  shown  that  the  three  curves  should 
share  the  critical  point.  For  polymer  pairs  having  plw  (and  p2w) 
different  from  1000  the  curves  remain  the  same  if  X  values  are 
suitably  rescaled.  This  can  be  seen  from  the  fact  that  the 
structure  of  the  Flory-Huggins  equation  remains  unaltered  when 
all  the  molecular  lengths  are  reduced  by  a  constant  factor  and  X, 
at  the  same  time,  is  multiplied  by  the  same  factor. 

The  most  interesting  feature  of  Fig.  1  is  that  as  the 
polydisperslty  Index  8  is  Increased,  the  cloud  point  curves 
become  flatter.  In  other  words,  at  compositions  away  from 
4>.|  s  0.5,  phase  separation  occurs  more  readily  (l.e.,  at  a  lower 
value  of  X)  as  the  polymers  become  more  polydisperse.  To 


understand  why  this  happens,  we  analyze  the  composition  of  the 
minority  phase  (or  phase  A)  which  just  begins  to  come  out  as  a 
separate  phase  at  the  cloud  point.  Fig.  2  shows  the  calculated 

A 

values  of  4>  -j  vs.  x  (the  so-called  shadow  curve)  for 
p1w  =  p2w  =  1000  and  =  %  =  2*°*  The  cl°ud  point  and  the 
spinodal  curves  are  also  shown  there.  The  horizontal  line  ties  a 
point  on  the  shadow  curve  with  a  point  on  the  cloud  point  curve, 
indicating  that  these  two  points  correspond  to  the  compositions 
of  the  coexisting  phases  at  the  given  value  of  x.  Fig.  3  shows 

A 

the  weight-average  p1w  of  polymer  1  in  phase  A  separated  out  from 
the  principal  phase  of  composition  Thus,  when  a  mixture 

containing,  say,  20%  of  polymer  1  is  brought  to  the  condition  of 
incipient  phase  separation,  the  minority  phase  that  is  being 
formed  contains  68.67%  of  polymer  1,  instead  of  80%  as  would  be 
expected  if  both  polymers  were  monodisperse.  The  polymer  1  in 
phase  A  consists  predominantly  of  higher  molecular  weight 
fractions  of  the  original  polymer  so  that  its  weight-average  DP 
is  now  equal  to  1853.  Thus,  this  fractionation,  in  favor  of  the 
higher  molecular  weight  end  of  the  distribution,  is  responsible 
for  the  separation  occurring  prematurely  at  a  value  of  X  smaller 
than  otherwise. 

For  the  polymers  having  the  Schulz-Flory  distribution 
considered  here,  the  polydlsperslty  index  of  the  polymers 
separating  out  in  the  minority  phase  is  not  different  from  that 
in  the  principal  phase.  This  is  so,  since  from  equations  (4), 
(7)  and  (12),  we  have: 

wA(p)  =  (ak+1pk/k!)  exp  t-(a-  o)p)  (17) 
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showing  that  the  index  g  =  1  +  1/k  is  unaltered.  The  same 
conclusion  is  reached  even  when  one  of  the  components  involved  is 
a  low  molecular  weight  solvent.  Thus,  in  the  case  of  a  polymer 
having  the  Schulz-Flory  distribution,  the  width  of  the 
distribution  cannot  be  made  narrower  by  fractional 
precipitation. 

Figs.  4  and  5  show  the  cloud  point  curves  calculated  for  the 
cases  where  p1w  =  p£W  =  1000  but  8^  6 2 •  The  cloud  point  curves 

develop  assymmetry  as  the  polydispersity  indices  of  the  two 
polymers  deviate  from  each  other.  The  distortion  of  the  curve 
away  from  the  symmetrical  shape  is  in  the  same  direction  as  one 
would  have  found  if  one  had  made  the  pw  of  one  of  the  polymers 
larger,  instead  of  making  its  distribution  broader. 

One  of  the  most  practical  ways  of  determining  X  parameters 
between  polymer  pairs  is  to  determine  the  cloud  points 
experimentally  and  then  find  the  best  value  of  X  which  brings  the 
observed  values  in  agreement  with  the  cloud  point  curves 
calculated  in  accordance  with  the  Flory-Huggins  expression.  In 
such  a  procedure,  the  effect  of  polydispersity  is  normally 
neglected  because  of  computational  difficulties.  The  present 
results  show  that  unless  the  polymers  employed  have  fairly  narrow 
distributions,  the  X  values  thus  obtained  could  entail  an 
appreciable  error.  When  the  two  polymers  are  symmetrical  with 
respect  to  pw  and  B,  the  error  is  more  likely  to  affect  the 
temperature  coefficient  of  the  X  parameter,  while  when  the  two 
polymers  differ  in  their  polydispersity ,  the  error  is  more  likely 
to  affect  the  concentration  dependence  of  the  X  parameter. 
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In  the  analysis  of  the  cloud  point  curves  of  polymer-solvent 
systems,  it  was  shown^  that  a  computational  difficulty  arises 
when  the  molecular  weight  distribution  is  a  "divergent"  type  (for 
example,  log-normal  distribution),  and  also  that  in  the  vicinity 
of  the  critical  point  a  more  complicated  phase  relations, 
sometimes  leading  to  a  coexistence  of  three  phases,  could  be 
encountered.  Analogous  situations  are  expected  to  arise  in  the 
case  of  polymer-polymer  systems  as  well,  but  we  have  not 
investigated  these  subtleties  in  this  work. 
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Figure  1.  Cloud  point  curves  calculated  for  the  symmetrical 


Figure  2. 


Figure  3. 


Figure  4. 


Figure  5. 


cases:  p1w  =  p2w  =  1000  and  81  =  62  =  1.0,  1.1  and 

2.0  as  indicated. 

The  cloud  point  curve  (labeled  £  and  the  spinodal 
curve  (broken  line),  calculated  for  the  case  where 
Plw  =  P2W  =  1000  and  ^  =  £2  =  2.0.  The  curve  labeled 

A 

<P  ■)  is  the  "shadow"  curve  giving  the  composition  of  the 
incipient  minority  phase  A.  The  horizontal  lines  tie 
pairs  of  points  giving  the  composition  £1  of  the 
principal  phase  and  the  composition  £^  of  the  minority 
phase  which  coexist  at  the  cloud  point. 

The  ordinate  gives  the  weight-average  DP,  P^‘K,  of 
polymer  1  in  the  incipient  minority  phase  A  which 
coexists  with  the  principal  phase  of  composition  £^, 
calculated  for  the  case  where  p1w  =  p2w  =  1000  and 
B1  :  B2  s  2.0. 

Cloud  point  curves  for  p1w  =  p2w  =  1000,  ^  =  2.0  and 

B2  =  1.0,  1.1  and  2.0  respectively.  As  the 

polydispersities  of  the  two  polymers  deviate  more  from 
each  other,  the  curves  become  more  unsymmetrical .  It 
is  as  if  the  molecular  weight  of  the  broader  polymer 
has  become  much  higher,  even  though  in  fact  the 
weight-average  DP's  of  the  two  polymers  are  the  same. 
Cloud  point  curves  for  plw  =  p2w  =  1000,  61  =  1  and 

$2  -  1*0,  1.1  and  2.0  respectively.  See  the  caption 
for  Figure  !». 
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